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To do so a pion potential is extracted from the in-medium dispersion rela- 
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^ ' tion and used in QMD calculations and thus we take care of both, real and 
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imaginary part of the pion optical potential. The potentials are determined 



■ from different sources, i.e. from the A-hole model and from phenomenologi- 

■ cal approaches. Depending on the strength of the potential a reduction of the 

■ anti-correlation of pion and nucleon flow in non-central collisions is observed 
' as well as an enhancement of the high energetic yield in transverse pion spec- 

tra. A comparison to experiments, in particular to pf-spectra for the reaction 

^ ■ Ca+Ca at 1 GeV/nucleon and the pion in-plane flow in Ne+Pb collisions at 

^ ■ 800 MeV/nucleon, generally favours a weak potential 
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I. INTRODUCTION 

The main motivation to study heavy-ion reactions at intermediate energies is to extract 
some information on nuclear matter under extreme conditions, i.e. at high densities and/or 
high temperatures. Besides nucleonic observables like the rapidity distribution and flow 
observables as, e.g., the bounce-off and squeeze-out also mesons |T|,0] (vr, K~^, i] etc.) emitted 
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from the reaction zone can be probes of excited nuclear matter. Pion yields and spectra 
turned out not to be sensitive to properties of the nucleon-nucleon interaction and 
difficult to interpret. Due to their strong interaction with the nuclear environment pionic 
observables at the freeze out are the result of complex creation and rescattering processes. 
However, pions can provide information about the resonance production in compressed and 
excited nuclear matter, in particular with respect to the A(1232) resonance which is the 
dominant production channel. In high energetic collisions with heavy nuclei the creation of 
"resonance matter", i.e. a state with a density of A-resonances comparable to the nuclear 
saturation density has, e.g., been discussed in this context ||l]J^. 

Whereas total pion yields and low energetic p^spectra can be well explained by present 
theoretical transport approaches the subthreshold production of high energetic pions is 
a question of current debate. Such spectra have, e.g., been measured by the TAPS and 
the KaoS collaborations and at the Fragment Separator at GSI. Various processes 
can contribute to the creation of subthreshold pions. One is the accumulation of energy 
by multistep collision processes, e.g. AA-scattering. But also heavier resonances as the 
A^*(1440) can be excited and decay into high energy pions. The analysis in particular of 
these high energy p^-spectra is supposed to yield information about the role of resonances in 
heavy ion collisions p. However, both cases are rare and require a sufficiently high density of 
excited resonances, a state which mainly occures in heavy systems. In lighter systems which 
have also been measured (Ar+Ca by TAPS and Ni+Ni in Ref. @]) multistep processes 
and the creation of higher resonances are less probable and it is not clear if these phenomena 
are sufficient to explain the experimental data. 



However, most theoretical approaches P,10| included the interaction of the pions with the 



surrounding nuclear medium only by collision processes, i.e. parametrizing the imaginary 
part of the pion optical potential. In the present work we include also the real part of the pion 
optical potential which influences the pion propagation through the nuclear medium. We 
thus take care of the full in-medium pion optical potential. We investigate the influence on 
Pt-spectra with pion momenta of several hundred MeV/c. This allows to test the in-medium 



dispersion relation in the high energy range. 

To investigate the low energy range of the dispersion relation we consider the momenta 
of the pions in the collective in-plane pion flow which have values of some 10 MeV/c. The 
pion flow as been measured, e.g., by the DIOGENE group for the system Ne+Pb at 800 
MeV/nucleon |TT|. The prominent observation was there an always positive pion flow in 
contrast to the typical S-shape of the nucleon flow. This fact is probably due to the projectile- 
target geometry in this highly asymmetric system and the reabsorption of pions by the target 
matter, i.e. the so called shadowing effect. In symmetric systems this shadowing effect 
may even lead to a complete anti-correlation of pion and nucleon flow [p|,p^. Although 
conventional calculations, i.e. without a pion potential, are able to qualitatively reproduce 
the DIOGENE data [0,0 in particular this (anti) correlation of flow reacts sensitive on the 



in-medium effects taken into account in the present work. 

The knowledge of the pion optical potential from elastic pion-nucleus scattering is, how- 
ever, restricted to nuclear densities at and below saturation and relatively small energies 



In intermediate energy heavy ion reactions up to 2 GeV incident energy per nucleon baryon 
densities of three times saturation density are reached and transverse momenta of about 1 
GeV/c are detected in pion spectra. In this range the real part of the pion potential is almost 
unknown. Thus, one has to extrapolate the pion dispersion relation to the ranges relevant 
for heavy ion collisions. As done in Ref. [0 we apply two models, i.e. the A-hole model 



||T5| , p!6[l and a phenomenological ansatz suggested by Gale and Kapusta |T3- The application 
of these models allows to investigate the influence of the possible boundary soft 
(A-hole) and a rather strong phenomenological potential. Moreover, since the in-medium 
corrected dispersion relation suggested by Gale and Kapusta is close to the hypothetic pion 
condensation limit it is a question of interest whether the boundary conditions occuring in 
intermediate energy heavy ion collisions would allow such a scenario. In addition we apply a 
third potential which is a modification of the parametrization of Gale and Kapusta and lies 
in strength between the two former cases. In Ref. [0 we have demonstrated the influence 
on pionic flow observables in heavy ion collisions. Here we want to derive quantitative state- 



merits and to distinguish between the different approaches by comparison to experimental 
data, i.e. to p^spectra of the emitted pions and the pion in-plane transverse flow produced 
in non-central collisions. 

The paper is organized as follows: First we give an overview of the treatment of pions 
in the Quantum Molecular Dynamics (QMD) approach, in particular with respect to the 
resonances included and the collision processes. Next we discuss the dispersion relation and 
the resulting pion potentials and then compare the results to experiment. Finally a summary 
and an outlook is given. 

II. PIONS IN QUANTUM MOLECULAR DYNAMICS 

QMD is a semiclassical transport model which accounts for relevant quantum aspects like 
the Fermi motion of the nucleons, stochastic scattering processes including Pauli blocking 
in the final states, the creation and reabsorption of resonances and the particle production. 
A detailed description of the QMD approach can be found in Refs. [0,|T^. 

Each baryon is represented by a Gaussian wave packet with a fixed width L. The tem- 
poral evolution of the centroids of these wave packets is described by the classical equations 
of motion generated from an A^-body Hamiltonian 

d^jj _ OHb 9q£ _ OHb 
dt dcii ' dt dpi 

with 

hs = j:M^+It. {u.^+urj'-'+w') ■ (2) 

The Hamiltonian, Eq. (H), contains mutual two- (and three-) body potential interactions 
which are finally determined as classical expectation values from local Skyrme forces Uij 
supplemented by a phenomenological momentum dependence and an effective Coulomb in- 
teraction 
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U., = a(P^]+p(P^y + 5ln' (e|p, - p,f + l) ^ 
\PoJ \PoJ ^ ^ Po 



(3) 



^r=(T) ^7— ^^-/(^^^) (4) 



where pij is a two-body interaction density defined as 

(47rL) 2 

and er/ is the error function. The parameters a, j3, 7, 5, e of the Skyrme interaction, Eq. @, 
are determined in order to reproduce simultanously the saturation density (po = 0.16/m~^) 
and the binding energy {Eb=-^Q MeV) for normal nuclear matter for a given incompress- 
ibility as well as the correct momentum dependence of the real part of the nucleon-nucleus 
optical potential [|1^,^. In the present calculations we apply the standard parametrization 
of a hard/soft equation of state (K=380/200 MeV). The Yukawa-type potential Uj^^^ in Eq. 
(0) mainly serves to improve the surface properties and the stability of the initialized nuclei. 

Analogously to the baryons, Eq. also the mesons obey Hamilton's equations of 

motion. The Hamiltonian of the pions, i.e. the sum of the respective single particle energies 
tUj, is represented in a mean field form 



H. = T.^^ = T. ( Vp' + < + R^v7'ip^^ q^) ) (6) 

i i 

where N^^ is the actual number of pions. The dependence of the pionic mean field, i.e. 
the real part of the pion optical potential ReV°^'^, on the pion coordinates originates from 
the medium effects in dense nuclear matter. Such medium effects are naturally expressed 
through a density dependence and ReV^^ takes the form 

i?ey7*(p„ q.) = ReV:^' (p., ^^(q.)) (7) 

where pB is the respective baryon density. In conventionell approaches the Hamilto- 

nian, Eq. (|^), is taken as that of the vacuum, i.e. ReV°^*' is set equal zero. 

Hard core scattering of the particles is included by the simulation of the collision pro- 
cesses by standard Monte Carlo procedures. The collisions probabilities are determined by 



a geometrical minimal distance criterium d < ycrtot/'^ weighted by the Pauli blocking fac- 
tors of the final states [|1^,^]. For the inelastic nucleon-nucleon channels we include the 
A(1232) as well as the A^*(1440) resonance. In the intermediate energy range the resonance 
production is dominated by the A, however, the A^* yet gives non- negligible contributions 
to the high energetic pion yield The resonances as well as the pions originating from 
their decay are explicitely treated, i.e. in a non-perturbative way and all relevant channels 
are taken into account. In particular we include the resonance production and rescatter- 
ing by inelastic NN collisions, the one-pion decay of A and N* and the two-pion decay of 
the A^* and one-pion reabsorption processes. These are summarized in Table 1 where the 
isospin dependence of the various channels is suppressed for the sake of simphcity; which is, 
however, taken into account in the calculations. 

For the cross sections of the inelastic NN channels we adopt the parametrizations of Ref. 
which have been determined from one-boson-exchange amplitudes in Born approxima- 
tion. The lifetimes of the resonances are determined through their energy and momentum 
dependent decay widths 

(1 + a2|p|^)(a3 + IpI^) 

which originates from the p-wave representation of the resonances. In Eq. (P) p is the 
momentum of the created pion (in GeV/c) in the resonance rest frame. According to Ref. 



21 the values ai=22.83 (28.8), a2=39.7 and a3=0.04 (0.09) are used for the A (A^*) and 



the bare decay widths are taken as r^= 120 MeV and r^*= 200 MeV. The reabsorption 
cross sections (vrA^ A, N*) are determined assuming a Breit-Wigner distribution for the 



^ p2(M^r^)2 {s - MIY + {MnVaf ^ > 

where R stands for a A or A^* and p is the cm. momentum. In contrast to, e.g., Ref. 
]TI| the momentum dependent resonance width, Eq. (^ is used in Eq. (|^). The resonance 
rescattering cross section {N{A, N*) —>■ NN) is determined from detailed balance. To 
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take care of the proper phase space available for this channel we take the finite width of 



the resonances into account as proposed in Ref. In particular for pions near the A- 

threshold this procedure leads to a significant enhancement of the rescattering cross section 
and the low energetic pion yield is suppressed by about 30%. Thus we are able to reasonably 
reproduce the total pion multiplicities with results similar to Ref. [^. 

In the presence of a pion potential this has to be taken into account in the energy balance 
in order to ensure energy conservation in the collision processes, i.e. the relation 



Pl + M| = J (PR - p.)2 + Mj, + Jpl + ml + ReV:^' (p., ps) (10) 



has to be fulfilled. Concerning the pion reabsorption Eq. ( [T0| ) is exactly fulfilled; in the 
case of a resonance decay this procedure is more involved since the decay takes place in the 
resonance rest frame. The potential of the created pion depends, however, on its momentum 
in the rest frame of the colliding nuclei. Thus the pion momentum is determined iteratively 
until Eq. (0) is fulfilled with an accuracy better than 0.5 MeV which corresponds to an 
energy violation of less than about 0.04%. 



Although the elementary cross sections p2[ are parametrizations of the free cross sections 



medium corrections to the imaginary part of the pion optical potential are included firstly 
by the enhancement of the rescattering/reabsorption probability which is proportional to 
the nuclear density, i.e. the available scattering partners, and secondly by the Pauli blocking 
in the final states according to the respective phase space occupancy. Via Eq. (0) also the 
real part of the pion potential gives corrections to the imaginary part. Thus the mean free 
path of the pions is significantly reduced and a consistent treatment of real and imaginary 
part of the optical potential is achieved. Furthermore, the momentum dependence of the 
nuclear mean field, Eq. (^, results in a (non-relativistic) effective mass of the baryons. 

III. THE PION POTENTIAL 



In contrast to the imaginary part of the pion optical potential which is well known from 
inelastic pion-nucleus scattering the knowledge of the real part extracted from elastic pion- 
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nucleus scattering is rare and restricted to low energies [|T4|. Data are available only for 
momenta up to about 2m^c. However, in intermediate energy heavy ion reactions pion 
momenta of more than 1 GeV/c can occure and the high energy tails of the pion spectra 
are of particular interest since they are supposed to yield information about the resonance 
production in hot and dense nuclear matter and corresponding in-medium effects 
Furthermore, nuclear densities up to three times saturation density can be reached whereas 
the elastic pion-nucleus scattering is restricted to densities around saturation and below. 
Thus, one has to extrapolate the dispersion relation to the energy range relevant for heavy 
ion collisions. However, such an extrapolation is affected by large uncertainities. To cover 
the unknown ranges we apply two contrary approaches, i.e. a microscopic and a phenomeno- 
logical ansatz for the pion dispersion relation. 



The microscopic ansatz is based on the perturbation expansion of the A-hole model [jl5 
The summation of the A-hole polarization graphs results in a pion self energy H entering 
into the in-medium pion dispersion relation 

o^(p)2 = p2 + m2 +H(u;,P,Pb) . (11) 

Here the self energy depends on the energy uj and the momentum p of the pion and on the 
nuclear matter density ps- In the framework of the A-hole model the self energy is finally 
given in the form 

na;,p,pB =1 — n V 12 



/ \ 2 

8 / /a \ c^a(p)pb 



9 - ^^(p) 



with 

x(^,P,Pb) ■ 

UJA = \jMl + p2 - M;v • 

The parameters entering into Eq. (|1^), in particular the ttA^A coupling constant /a and 
the correlation parameter g' are taken in consistence with the OBE parameters of Ref. ||2^ 



and a consistent treatment of the real and imaginary part of the pion optical potential is 
achieved. 
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Although this approach works reasonably well at low baryon densities and low pion 
momenta |[T^ one has to be careful when appying Eq. (|1^) to heavy ion reactions. The 



self energy obtained from the A-hole model includes beside of excitation of AA^~^ states 
also short range correlations of these states. In this approximation one neglects, however, 
terms of higher order which are necessary to prevent the system to undergo a phase 



transition to the so called pion condensation. Furthermore, Eq. ( [Til) yields the wrong 
boundary conditions for high energy pions when they pass through the surface of nuclear 
matter into the vacuum, in particular a "quasi-pion" expressed by a A-hole state has to 
convert asymptotically to a real pion which can be detected. A common practice to avoid 
these unphysical features is to mix the two solutions of the dispersion relation, Eq. (pTj), 
i.e. the pion-like and the A-hole-like branch {^1,0,2) in order to obtain an effective pion 



dispersion relation 25 



UJeff = ZiVLi + Z2^2 , Zi + ^2 = 1 . (13) 

The probabilities for the quasi-pion to sit on the pion-like branch [Zi) or the A-hole-like 
branch (Z2) are determined from the condition 



I ^ Z^ ^ Z2 

^2 _ p2 _ ^2 _Y[ oj"^ — Vl{ uj"^ — VL2 



2 _ r,2 _ ^2 _ TT ,.,2 _ 02 + ,.,2 _ 02 ^^^^ 



and thus the physical boundary conditions are fulfilled. Replacing oj in Eq. (0) by cUg// one 
obtaines an effective pion self energy He//- The real part of the optical potential is obtained 
from the pion wave equation, i.e. the Klein-Gordon equation as 



ReVr{v.PB) = ^'''^'"''''^'^^^ . (15) 

The density and momentum dependence of the dispersion relation relation are shown in 
Fig. |I] and it is seen that such a construction leads to a strong softening of the in-medium 



effects compared, e.g., to the original dispersion relation ||2^. Consequently, the medium 
dependence of the resulting potential, Eq. (|1^), (called Pot.l in the following and shown in 
Fig. |I|) is moderate and the attraction of the potential is weak. Hence one can't be sure 
that the modified dispersion relation still represents the true pion-nucleon interaction. 



As already done in Ref. |T3| we also consider a phenomenological ansatz suggested by 



Gale and Kapusta |T^. The dispersion relation then reads 



uj{p) = J{\v\-Po? + ml-U (16) 



U = ^pl + ml- (17) 
mo = m^r + 6.5(1 — x^^)'m.„ (18) 
pI = {I - xYm\ + 2mom^(l - x) . (19) 

A phenomenological medium dependence is introduced via x = e""*^''^/''") with the parameter 
a = 0.154 and po the saturation density in nuclear matter (here po = 0.16/m~^). The form 
of Eqs. (|16| - [18|) is motivated by the following constraints : 



i) The group velocity duj/dp may not exceed the velocity of light. 

For high energetic pions many-body effects should be of minor importance and uj 
should in leading order be proportional to p for p — > oo. 



m 



IV 



Corresponding to a strong p-wave interaction the energy should have a minimum. 

Due to a weak s-wave interaction medium effects are weak for pions at rest with respect 
to the surrounding medium. 

The observation of pionic atoms implies that uj ~-20 MeV for ps = Po and p ^ 2111^^0. 
(vi) Pion condensation can only appear at infinite density. 
The dispersion relation as well as the corresponding potential (called Pot. 2 in the following) 



are shown in Fig. |I|. Here ReV^^ is extracted from Eq. (|16]) as ReV^^ = — y + m^. 
It turns out that the potential is much stronger and its medium dependence is much more 
pronounced than in the case of the modified A-hole dispersion relation. In this context we 
want to mention that this effect is to a large extent due to the mixing of the pionlike and the 
A-hole-like branches in Pot.l since a potential extacted according Eq. (|TB|) from the pure 
pionlike branch given in A-hole model is of the same magnitude as Pot. 2 ||26|. Thus, Pot.l 
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and Pot. 2 can be looked as the boundary cases of a rather weak and a very strong pion 
potential. Furthermore, the phenomenological dispersion relation 2 is closer to the limit 
where pion condensation can occur. 

In order to test the pion dispersion relation it is therefore a natural step to suggest a third 
parametrization (Pot. 3) which lies in between the former ones. To do so, we decrease the 
power x^^ in Eq. (|18D to x^. The constraintes (i)-(iv) remain unrendered by this modification 
except of point (v) where now a value of +10 MeV is obtained instead of -15 MeV in the 
former case. However, we do not consider this as a severe drawback since the application of 
the potential to heavy ion collisions should not be affected therefrom. In addition, also the 
A-hole model potential does not yield a bound quasi-pion for these values of pb and p. 



IV. RESULTS 

First we consider transverse momentum p^-spectra for the system ^°Ca+^°Ca at 1 GeV 
incident energy per nucleon. In particular we compare to 7r°-spectra recently measured by 
the TAPS collaboration 0. Here the vr^'s have the advantage that they are not distorted 
by Coulomb effects. The calculations are impact parameter averaged and a rapidity cut of 
—0.2 < yc.m. ^ 0.2 has been applied which takes into account the detector acceptance. The 



data shown in Figs. |^ and ^ have been obtained with this cut ||2^ and are slightly enhanced 
with respect to the results shown in Ref. [§] where a larger cut has been used. In fig. 2 we 
first investigate the influence of the nuclear equation of state and therefore compare a soft 
(SMDI) and a hard (HMDI) momentum dependent Skyrme force, see Eq. (^. It is seen that 
the dependence of the pion spectrum on the nuclear EOS is moderate and the agreement 
with the data is reasonable, however, not overwhelmingly good. We want to mention that 
the results obtained with the hard EOS are close to those of Ref. 0. As also found there 
the low pt-spectrum is slightly overestimated whereas the high pt range is underestimated 
by about a factor of 3. In Ref. the lack of high energy pions is even more pronounced 
than in the present calculations. This is probably due to the A^*(1440) resonance which was 
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not taken into account there. 

Next we turn to the influence of the in-medium pion potential. In Fig. ^ we show the 
7r°-spectra for the same reaction as in Fig. ^ now including the pion-nucleon interaction as 
discussed in the previous section. All calculations are performed with the soft momentum 
dependent Skyrme force. It is seen that the weak A-hole potential (Pot.l) has nearly no 
influence on the prspectrum which is in agreement with the results of Ref. ||2^ where a 
similar potential has been investigated. In this context we want to mention that concerning 
the low energy pions the present results (and those of Ref. |^) stand somehow in contra- 
diction to the BUU calculations of Refs. p5| , |26| since there a significantly underestimated 
low energy pt yields was obtained when the pion potential was neglected. However, in these 
works comparisons to other measurements have been performed [p8| , p9| . In our opinion 
the discrepancies are rather caused by the analysis of the data than by a slightly different 
treatment of the pion creation and annihilation processes outlined in section 2. 

In the case of the phenomenological potentials 2 and 3 the situation is completely different 
now. With Pot. 2 the low p^-range is slightly enhanced, however, the high energy pions are 
overestimated by nearly one order of magnitude. Such an enhancement of low energy pions 
has been observed in Ref. [03] where an effective pion potential has also be determined 



within the A-hole model which is, however, apparently stronger than Pot.l in the present 
work. This complex behavior can be understood by the strong attraction of the respective 
potential which forces the pions to follow the trajectories of the nucleons. Most pions get 
bound by the stopped participant matter resulting in an enhancement of the low pt-jield. 
On the other handside, pions which are bound by the spectator matter are driven out 
to high transverse momenta by the nucleonic flow and thus the high pt-range is strongly 
enhanced. This effect is diminished when the potential is weaker (Pot. 3). Here the very 
high energy yield {pt > 0.7GeV/c) is still slightly overestimated but also the low pf-range is 
now overestimated by about a factor of 3. The latter is undestandable since for low momenta 
{p < l.Sm^c) Pot. 3 is even stronger than Pot. 2, see zoomed region in Fig. |T]. 

The above observations are clearly reflected in the pion transverse flow. In Fig. ^ 
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we compare the pion and the baryon in-plane transverse flow for the reaction Ca+Ca at 
1 GeV/nucleon. The results are impact parameter averaged and scaled to the pion and 
nucleon mass, respectively. By this scaling pion and nucleon flow are of the same order 
of magnitude. It is seen that without any pion-nucleon interaction the flow of pions and 
nucleons is clearly anti-correlated. The anti-correlation originates from the absorption of 
the pions by the participant matter which produces the nucleonic flow. This is known as 
the shadowing effect. In the presence of the weak A-hole potential (Pot.l) the pion flow 
remains nearly the same. In the case of Pot. 2 the pion flow switches from an anti- to a 
correlation with the nucleon flow, i.e. the pions are forced to follows the trajectories of the 
spectator nucleons by the strong attractive potential. Further it is seen that the strength 
of Pot. 3 has just a magnitude where the attraction and the shadowing effect are nearly 
completely counterbalanced and the resulting flow is around zero. This overall behavior is 
in good agreement with the interpretation of prspectra given above. 

In summary none of the cases under consideration is able to reproduce the experimental 
spectrum over the entire range of energy with a satisfying high accuracy. Both, the low 
as well as the high p^-range react sensitive on the pion potential. Pot. 2 seems to be too 
attractive and can be ruled out by the comparison to the data. However, the present results 
indicate that the inclusion of higher resonances may be not sufficient to explain the role of 
high energy pions but furthergoing medium effects should be taken into account. 

Next we investigate the collective in-plane transverse flow of pions which has been mea- 



sured for the system Ne+Pb at 800 MeV/nucleon by the DIOGENE collaboration |TT|. There 
the most striking results was the observation of a positive pion flow also for backwards ra- 
pidities, i.e. the pion flow is partially anti-correlated to the nucleon flow. In particular here 
we consider a semi-central reaction at impact parameter b=3 fm. For a comparison with 



the data we included the experimental detector fllter cuts given in Ref . |Tl[ . In addition, we 
simulated the reconstruction of the estimated reaction plane as it was done in the experi- 
mental analysis. In contrast to the theoretical calculation where the true reaction plane is 
known a priori in the experiment the transverse in-plane vector is estimated for every event 
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by 

Q = E^.P^i (20) 

3 

where the sum runs over all detected protons weighted by their relative rapidity ujj = 
Hi— < y > with < y > the mean rapidity of the total system. The estimated in-plan pion 
momentum = ■ Q is then obtained by the projection on the unit vector Q. 

Fig. ^ shows the corresponding vr^ rapidity distribution obtained with the various pion 
potentials and including the detector cuts. The calculation without pion potential is thereby 
in a good agreement with the result of Ref. |[10|. The influence of the pion potentials 1 and 
2 on the longitudinal flow is relatively weak, only in the case of Pot. 2 a slight enhancement 
of the backward, i.e. the target-like, rapidity distribution is observed. Applying Pot. 3 we 
observ a strong suppression of the detected 7r+ yield over the entire rapidity range. This 
effect is also reflected in the total yields given in Table 2; is, however, not so pronounced since 
the total yield is only suppressed by about 10%. In generel the usage of a potential leads 
to a slight reduction of the total yields which was also found in Refs. P3| , |2B|] . Concerning 
Pot. 3 the effect seen in Fig. |^ seems to be somehow an artefact of the detector geometry. 
Nevertheless, the dynamics of the pions is significantly changed and slow pions near the 
A-threshold get in particular strongly bound and and are preferrentially reabsorbed by the 
spectator matter. 

In Fig. |] the corresponding in-plane flow per pion is shown in units of the pion mass. 
Here the procedure to reconstruct the experimental reaction plane has been performed. It 
turns out that in the case of a vanishing or a weak potential (Pot.l) the flow is always 
positive and in a fairly good agreement with the data for values around mid-rapidity. In the 
high forward rapidity range we, however, underestimate the data by about a factor of two. 
Since the results are strongly distorted by the large asymmetry of the considered system 
and the reconstruction of the reaction plane this effect can be due to a lack of sufficiently 
good statistics in this range. The positive flow at backward rapidities is explained by the 
reabsorption of the pions by the large target and thus is a consequence of the shadowing 
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effect. In the case of the strong Pot. 2 the shadowing effect is overcompensated by the 
attraction of the potential and the flow shows a definit change of sign from negative to 
positive values at midrapidity. Similar as in Fig. ^Pot.3 shows the same behavior which is, 
however, not so pronounced. Both, Pot. 2 and 3 yield significantly too less positive flow over 
the entire rapidity range. 

These observations are also reflected in the ratio R between the numbers of pions with 
positive and negative values oipx displayed in Table. 2. Since we did not explicitely include 
the Coulomb interaction in the propagation of the pions, Eqs. (|l]) and (^, the different 
charge states are only taken into account via the isospin dependence of the respective cre- 
ation/absorption channels. Thus, for a fair comparison to the measured ratio a mean value 
(corrected for the respective total yields) of vr"'' {R = 1.42) and vr^ (i? = 1.30) which is about 
R = 1.34 should be compared. Then the theoretical positive abundancy is a reasonable 
agreement with the experiment for the calculations without pot. and with Pot.l. In the 
case of Pot. 2 and Pot. 3, however, R is close to unity. 

V. SUMMARY AND CONCLUSIONS 

We investigated the influence of medium corrections to the pion dispersion relation on 
the pion dynamics in heavy ion collisions. This was done by the introduction of a pion- 
nucleus potential through which the pions propagate between their collisions and refers 
to the real part of the pion optical potential. The imaginary part is included by the non- 
perturbative description of elementary rescattering and reabsorption processes and is directly 
medium corrected via the respective densities of nucleons and resonances. Hence a consistent 
treatment of the real and imaginary part of the optical potential was achieved. 

The pion potential was extracted from the dispersion relation given in the A-hole model 
thereby mixing the pion-like and the A-hole-like branches in order to avoid some unphysical 
features of the model. This, however, leads to a strong softening of the medium effects. 
Thus we also applied a phenomenological dispersion relation with a pronounced medium 
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dependence. 

Studying the influence of such effects the different approaches have been subjected to a 
comparison to the experiment for two contrary observables, i.e. spectra and in-plane flow. 
Thereby one has to keep in mind that pions at the freeze out are remanents of a complex 
collision history, i.e. of multiple creation and reabsorption processes. The analysis of the 
DIOGENE flow data implies a preferrential abroption of pions by the spectator matter 
known as the shadowing effect which may even lead to an anti-correlation of the collective 
pion and nucleon in-plane flow. However, a strong attractive pion potential leads to a 
bending of the pions by the nucleons and thus was found to favour a correlation of the 
respective flow. The magnitude of this effect is directly proportional to the strength of the 
potential. Since the pions get bound by the participant matter (slow pions) as well as by 
the spectator matter (fast pions) this process is complex, i.e. an enhancement of the low 
Pt-range as well as of the high p^-range was observed in transverse momentum spectra. 

A comparison to experiments, i.e. the TAPS spectra and to the DIOGENE flow data, 
seems to rule out a too attractive optical potential. Hence the pions created in intermediate 
energy heavy ion collisions are far from the limit of pion condensation. The weak effective 
A-hole potential has nearly no influence on pionic observables as well as the dependence 
on the nuclear equation of state was found to be rather weak. However, in particular the 
spectra of high energetic pions react sensitiv on the dynamics and the present results would 
favour a pion potential which's attraction is weak at low energies but becomes more pro- 
nounced with increasing energy. In our opinion a furthergoing analysis of high energetic 
subthreshold spectra as, e.g., measured in Ref. |^ may help to clarify the question if conven- 
tional approaches as, e.g., energy accumulating multi-scattering processes or the creation 
of higher resonaces are sufficient to explain such data. The present results indicate that 
the inclusion of in-medium corrections to the real part of the pion optical potential is of 
essential importance for a correct description of pion djTiamics. In addition the analysis of 
the pion in-plane flow also in symmetric systems as, e.g. recently measured by the FOPI 
collaboration, will help to learn something about the pion dispersion relation from heavy 
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ion 



collisions. 
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FIG. 1. Upper Figure: The free pion dispersion relation (solid) is compared to various 
in-medium dispersion relations at nuclear densities pB = Po (respective upper curves) and pB = 2po 
(respective lower curves). The in- medium dispersion relation is extracted from the A-hole model 
(Pot.l, long-dashed) and the phenomenological parametrization of Ref. ||T^ (Pot. 2, dotted) as well 
a softer parametrization (Pot. 3, dashed) are shown. Lower Figure: Real part of the pion optical 
potential at densities pb = Po (upper curves) and pB = 2po (lower curves) extracted from the 
corresponding dispersion relations. The inserted figure zoomes the low momentum range. 
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FIG. 2. Influence of the nuclear EOS on the vr'^ p^-spectrum for the reaction '^^Ca+'^^Ca at 1 
GeV/nucleon. A hard (HMDI) and a soft (SMDI) momentum dependent Skyrme force have been 
used. The experimental data heave been measured by the TAPS Collaboration for the system 
40^r+40(^a at 1 GeV/nucleon |]. 
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FIG. 3. Influence of the in-medium pion potential on tlie vr'^ p^-spectrum for the same reaction 
as in fig. 2. The calculations are performed without (solid) and including a pion potential. The re- 
spective potentials have been obtained from the A-hole model (Pot.l) and by the phenomenological 



ansatz of Ref. |17] (Pot.2). Pot.3 corresponds to a modification of Pot.2. 
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FIG. 4. Influence of the in-medium pion potential on the in-plane tansverse pion flow as a 
function of the center-of-mass rapidity for the same reaction as in fig. 2. The calculations are 
performed without (solid) and with inclusion of a pion potential. The respective potentials have 
been obtained from the A-hole model (Pot.l, long-dashed) and by the phenomenological ansatz of 
Ref. ||T^ (Pot. 2, dotted). Pot. 3 (dashed) corresponds to a modification of Pot. 2. In addition the 
nucleon flow is shown (solid with diamonds). The results are scaled to the pion and nucleon mass 
and thus the transverse flow is given in units of rriT^c and mj^c, respectively. 
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FIG. 5. vr^ rapidity distribution for a semicentral (b=3 fm) Ne on Pb reaction at 800 
MeV/nucleon after applying the experimental detector filter cut. The calculations are performed 
without (solid) and including a pion potential. The respective potentials have been obtained from 
the A-hole model (Pot.l, long-dashed) and by the phenomenological ansatz of Ref. (Pot.2, 
dotted). Pot.3 corresponds to a modification of Pot.2. 
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FIG. 6. Transverse in-plane flow per pion sacled to the pion mass for a semicentral (b=3 fm) 
Ne on Pb reaction at 800 MeV/nucleon after applying the experimental detector filter cut. The 
calculations are performed without (solid) and including a pion potential. The same potentials as 
in fig. 5 have been used, i.e. Pot.l (long-dashed), Pot. 2 (dotted) and Pot. 3 (dashed). The data 
are taken from Ref. [pT|]. 
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NN ^ NA Ntt 

NN ^ NN* N* ^ Ntt 

NN ^ AA N* ^ Avr 

ATA ^ AA N* ^ Ntttt 

TABLE I. Inelastic scattering processes which are included in the present calculations. 



No Pot. Pot.l Pot. 2 Pot. 3 Exp. 





1.89 


1.86 


1.80 


1.59 




R 


1.26 


1.25 


1.03 


1.12 


1.34 



TABLE II. Total vr^ yield (without experimental filter) and the ratio R of detected tt^ and vr^ 
pions with positive to negative values of px for the reaction Ne on Pb at 800 MeV/nucleon and 
impact parameter b=3 fm. The experimental value for R is an average of the respective values for 



IT and vr^ given in Ref. |11|. 
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